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SIMULATION  STUDIES  OF  PARTICLE  ACCELERATION  POWERED  BY 
MODULATED  INTENSE  RELATIVISTIC  ELECTRON  BEAMS 


I .  Introduction 

Future  progress  in  accelerators  and  their  applications  may  depend 

critically  on  the  development  of  physical  mechanisms  capable  of  generating 

high  voltage  gradients.  It  has  been  shown  that  a  high  electric  field  can 

be  established  in  rf  structures  by  modulated  intense  relativistic  electron 

9 

beams  (MIREBs)  of  power  erreater  than  10  V  and  that  such  a  beam  may  be  used 

as  an  rf  source  to  power  an  accelerator,  obtaining  voltage  gradients  as 

1  2 

high  as  100  MeV/m  or  greater.  ’  Such  accelerators,  in  which  a  low  power, 

high  current  beam  interacts  via  a  metallic  structure  with  a  low  current 

beam  to  obtain  very  high  energies  have  been  suggested  by  an  number  of 

3 

authors,  including  Voss  and  Veiland,  in  addition  to  the  accelerator 
outlined  in  Ref.  2.  Vakefield  acceleration  has  recently  been  observed  in 

4 

experiments  carried  out  by  Figueroa  et  al. 

2  5 

Theoretical  discussions  of  accelerators  powered  by  MIREBs  ’  have 
suggested  unusual  properties  that  may  be  present  in  the  following  devices: 
Firstly,  the  demonstrated  conversion  of  the  high  dc  power  of  an  intense 
relativistic  electron  beam  (IREB)  to  high  rf  power  in  the  MIREB  by  the  use 
of  tuned  radial  cavities  implies  that  the  MIREB  may  be  coupled  to  an  rf 
structure  so  as  to  drain  significant  power  (>  1  GW)  from  the  beam  at  high 
efficiency  and,  secondly,  geometrical  effects  may  allow  for  sizeable 
variations  in  efficiency,  field  gradient,  and  coupling  between  the  high 
power  MIREB  and  the  rf  structure  with  small  changes  in  the  experimental 
parameters . 

In  the  present  paper  we  study  these  issues  via  an  axisymmetric 
particle  simulation  using  the  CONDOR^  code,  which  has  been  previously  and 
successfully  applied  to  the  physics  of  such  intense  beams. ^  The 
accelerator  configuration  to  be  studied  is  similar  to  that  outlined  in  Ref. 
2  and  is  pictured  in  Fig.  1. 

Manuscript  approved  December  28.  1988 
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(1)  An  annular  IREB  generator  injects  a  beam  of  radius  r^  =6.3  cm, 

current  I  =16  k_A,  energy  E.  .  =  500  keV  and  duration  T  =  150  ns  into  a 
o  inj 

drift  tube  of  radius  ry  =  6.8  cm.  The  IREB  is  guided  by  an  axial  field, 

B  =  10  kG. 
o 

(2)  The  IREB  is  fully  modulated  at  f  =  1.3  GHz  by  a  pair  of  tuned 
radial  cavities,  the  first  of  which  is  externally  driven  by  a  low  level  rf 
source  (magnetron).  The  modulation  region  is  immersed  in  the  axial 
magnetic  field. 

(3)  The  MIREB  is  guided  into  a  cylindrical  cavity  of  radius  9.6  cm. 

The  cavity  is  loaded  with  thin  disks  of  radius  9.0  cm  and  separation 
1.88  cm.  The  MIREB,  which  has  a  frequency  of  modulation  corresponding  to 
the  desired  mode  of  the  rf  structure,  is  terminated  at  the  first  disk.  A 
resonant  interaction  occurs  at  the  gap  defined  by  the  end  of  the  drift  tube 
and  the  first  disk  of  the  rf  structure,  transferring  energy  from  the  beam 
to  the  rf  structure. 

(4)  An  emitter,  located  on-axis  on  the  surface  of  the  first  disk, 

emits  electrons  when  the  fields  within  the  structure  reach  a  sufficiently 

high  value.  This  secondary  beam  is  then  accelerated  by  the  rf  fields, 

guided  by  the  axial  magnetic  field. 

The  modulation  stage  of  this  device  has  been  studied  in  some  detail 

5  7  8 

for  a  1.9  cm  radius  annular  beam  in  a  2.4  cm  radius  drift  tube  ’  ’  and  has 
been  successfully  repeated  at  I  =  16  kA  and  r^  =  6.3  cm  in  a  drift  tube  of 

9 

radius  6.8  cm. 

In  the  present  paper,  we  will  investigate  the  coupling  between  the 
modulated  beam  and  the  rf  structure  and  the  subsequent  acceleration  of  the 
secondary  beam  and  shall  proceed  as  follows.  In  Sec.  II,  we  give 
theoretical  background  on  the  expected  field  gradients  in  the  rf  structure, 
define  a  transformer  ratio  for  this  acceleration  scheme,  and  present 
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numerical  results  from  the  Superfish^  code  on  the  modes  of  the  rf 
structure.  In  Sec.  Ill,  which  contains  the  main  results  of  this  paper,  we 
will  simulate  particle  acceleration  and  will  see  that  power  in  excess  of  1 
GU  may  he  transferred  between  the  primary  and  secondary  beams.  Here  we 
will  consider  the  effect  of  geometrical  variations  on  the  beam-rf  structure 
coupling  and  on  the  transformer  ratio.  Section  IV  will  contain  a  detailed 
discussion  of  the  numerical  issues  that  effect  the  ability  of  these 
simulations  to  correctly  predict  experimental  results.  Section  V 
concludes . 


II .  Fields  in  the  RF  Structure 

The  process  of  energy  transfer  between  the  primary  and  secondary  beams 
in  this  accelerator  resembles  that  of  the  Wakefield  schemes  described  in 
Refs.  3  and  4  in  the  use  of  fields  excited  by  the  primary  beam  in  a  disk 
loaded  structure.  In  these  schemes,  the  two  beams  travel  colinearly  such 
that  the  transformer  ratio  is  defined  as  R  =  ^2^1  w^ere  Ei  is  the 
magnitude  of  the  decelerating  field  experienced  by  the  primary  beam  and 
is  the  accelerating  gradient  experienced  by  secondary  beam.  In  the  present 
case,  however,  the  interaction  of  the  primary  beam  with  the  rf  structure 
takes  place  only  as  the  beam  traverses  the  gap  near  the  first  disk  of  the 
rf  structure,  where  the  beam  is  terminated,  while  the  secondary  beam  is 
accelerated  along  the  entire  length  of  the  rf  structure.  The  transformer 
ratio  is  then  defined  as 


R 


<E  .  >L 
axis 

E  d 
gap 


(1) 
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where  d  is  the  gap  length,  L  is  length  of  the  rf  structure,  E^  is  the 

decelerating  field  in  the  gap  (assumed  spatially  constant)  and  <E  .  >  is 

SIX  1  s 

the  average  field  experienced  by  the  accelerated  secondary  beam  particles. 
This  geometry  has  been  modelled  as  an  interaction  between  a 

sinusoidally  varying  current  source  and  a  transmission  line  consisting  of  a 

2 

series  of  R-L-C  circuit  elements.  This  model  exhibited  many  features  that 
have  been  found  in  the  numerical  simulations,  but  such  a  model  has  a 
limited  predictive  capability. 

Some  insight  into  this  problem  may  be  obtained  by  assuming  that  the 

disk  structure  will  behave  like  a  resonant  cavity.  The  normal  modes  of 

this  cavity  may  be  solved  for  by  neglecting  the  interaction  region  at  the 
2 

gap.  For  the  purposes  of  this  discussion,  however,  we  may  consider  the 
only  fundamental  mode  of  a  disk-loaded  cavity  of  length  L  =  nX/2  where  X  is 
the  wavelength  of  the  rf  and  n  is  a  positive  integer.  In  this  case,  the 
z-component  of  the  electric  field  of  the  fundamental  mode  varies 
sinusoidally  along  the  axis  and  radially  as 

Ez(r)/Ez(r=0)  =  J0(kr)/JQ(r=0)  ,  (2) 

where  Jq  is  a  Bessel  function  and  k  =  2it/X.  In  the  analysis  of  Ref.  2,  it 
was  conjectured  via  a  heuristic  argument  that  the  ratio  of  the  field 
experienced  at  the  gap  by  the  primary  beam,  r  =  r^,  to  the  peak  field 
on-axis,  r  =  0,  would  follow  this  radial  variation.  This  suggests  that  the 
radial  position  of  the  primary  beam  in  relation  to  the  mode  structure 
within  the  rf  cavity  is  of  some  importance  for  the  strength  of  the 
interaction,  the  efficiency,  and  the  obtainable  transformer  ratio. 

The  normal  modes  for  a  given  axisymmetric  cavity  may  be  calculated 
numerically  by  using  the  Super  fish^  code.  The  Superfish  result  for  one 
such  cavity  is  shown  in  Fig.  2.  Here,  a  disk  loaded  cavity  of  length  L  =  X 
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is  used  and  the  gap  region  is  included.  Except  for  the  metallic  boundary 
condition  imposed  at  the  right-hand  vail,  this  geometry  closely  resembles 
that  of  Fig.  1,  where  the  right-hand  boundary  is  an  open  drift  tube  for 
which  the  1.33  GHz  cavity  mode  is  below  cutoff.  This  result  and  a  series 
of  similar  results,  where  the  location  of  the  right-hand  wall  was  varied, 
show  that  the  expected  cavity  mode  is  obtained. 


Ill .  Numerical  Simulations 

The  simulation  geometry  (Fig.  3)  consists  of  a  short  drift  tube  region 
with  radius  r  =  6.8  cm,  a  gap  of  length  d  =  1.57  cm  and  a  disk-loaded 
structure  of  length  L  =  22.2  cm  =  X,  where  X  =  c/f,  and  f  =  1.27  GHz  is  the 
frequency  of  the  accelerating  mode  of  the  cavity,  and  was  determined 
numerically. 

The  primary  beam  is  injected  from  the  left-hand  wall  with  radius 
rb  =  6. A  cm,  energy  =  2.0  MeV  and  current  Iin^(kA)  =  g(t) 

[16  +  8sin(2nf t) ] ,  where  g(t)  is  an  envelope  function  that  increases 
linearly  from  zero  to  unity  during  the  time  0  <  t  <  15  ns  and  remains 
constant  thereafter.  At  a  selected  time,  t  >  15  ns,  the  secondary  beam 
with  I 2  =  10  Amperes  and  E£  =  0.1  MeV  is  injected  continuously  from  the 
center  of  the  first  disk  and  is  accelerated  along  the  axis  by  the  rf 
fields.  Each  simulation  continues  until  t  =  30  ns. 

Note  that  in  order  for  the  cavity-mode  approximation  of  Sec.  II  to  be 


of  use  the  parameters  L,  v  and  T  must  be  such  that  L/v  <<  T,  where  v  is 

g  g  g 

the  group  velocity  of  E-M  radiation  within  the  disk-loaded  structure  and  T 


is  the  duration  of  the  primary  beam  pulse.  If  this  relation  is  not 


satisfied,  the  disk  structure  will  behave,  not  like  a  cavity,  but  like  a 


travelling  wave  tube.  In  the  numerical  geometry  of  Fig.  3,  we  have 
arranged  the  separation  between  the  outer  disk  edges  and  the  cavity  wall  so 
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that  v  /c  =  1.  In  this  case  the  condition,  L/v  «  T,  is  easily  satisfied 
g  g 

within  the  30  ns  duration  of  the  simulations. 

Several  differences  between  this  configuration  and  that  of  a  practical 
experiment  must  be  noted. 

(1)  In  a  practical  experiment  the  disk  structure  would  be  longer  so  as 
to  obtain  higher  energies  in  the  secondary  beam.  Another  difficulty  is 
that  the  high  group  velocity  of  the  E-M  waves  in  the  simulation  structure 
and  its  short  length  and  small  volume  would  allow  rf  fields  to  build  up  so 
quickly  that  they  might  reflect  the  primary  beam  in  an  actual  device. 

(2)  An  actual  device  would  have  support  rods  to  hold  the  disks  in 
place.  These  would  also  provide  a  path  for  the  dc  current  of  the  primary 
beam.  Because  such  supports  cannot  be  modelled  axisymmetrically  and 
because  we  require  a  dc  current  path,  we  inserted  a  center  conductor  to 
serve  this  purpose.  It  will  be  seen  in  Sec.  IV,  below,  that  the  presence 
cf  this  renter  conductor  docs  not  significantly  effect  the  results. 

(3)  In  the  simulation  geometry,  the  left-hand  boundary  is  a  metallic 
wall.  In  an  actual  device  and  in  Fig.  1,  this  boundary  is  an  open  drift 
tube,  for  which  the  1.27  GH^7  frequency  of  the  rf  field  is  below  cutoff.  The 
metal  boundary  of  the  simulations  will  have  the  similar  effect  of 
reflecting  incident  radiation  at  this  frequency,  but  is  clearly  not  the 
same . 

(4)  The  simulation  structures  are  defined  on  a  grid  such  that  the 
effective  skin  depth  of  the  material  is  one  grid  cell  (Ar  =  0.2  cm, 

Az  =  0.3133  cm),  making  the  cavity  extremely  lossy,  with  Q  of  order  10. 

A  typical  value  for  a  metallic  structure  is  of  order  1000. 

Figures  4  and  5  show  the  z-component  of  the  electric  field  plotted  vs. 
time  in  the  gap  and  on-axis,  respectively,  for  a  simulation  with  parameters 
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described  above.  The  plot  on-axis  is  taken  at  the  spatial  location  of  the 
peak  electric  field.  Ve  see  that  the  fields  increase  continuously, 
reaching  values  of  56.3  MV/m  at  the  gap  and  94.2  MV/m  on-axis  before  the 
simulation  is  halted.  The  plot  of  the  gap  electric  field  shows  evidence  of 
a  weak,  lover  frequency  mode  which  may  have  been  excited  by  the  increase  in 
dc  current  from  t  =  0  to  t  =  15  ns.  The  rf  cavity  mode,  as  expected,  is  a 
standing  wave,  varying  sinusoidally  in  z  and  as  a  Bessel  function,  Jp(kr), 
radially.  This  is  seen  in  Figs.  6  and  7  which  show  Ez  vs.  z  and  vs.  r, 
respectively,  at  fixed  time. 

For  the  simulation  shown,  the  secondary  beam  was  injected  continuously 

for  t  >  17  ns  with  =  10  A  and  E ^  =0.1  MeV  and  was  bunched  and 

accelerated  by  the  rf  fields.  This  acceleration  may  be  observed  in  Fig.  8, 

which  plots  particle  positions  in  phase  space,  ygc  vs.  z,  where  0  is  the 

-2  -1/2 

axial  particle  velocity  normalized  to  c  and  y  =  (1  -  0  )  .  The 

particle  positions,  plotted  at  fixed  time  at  intervals  of  0.2  ns,  show  a 
maximum  energy  increase  of  8.60  MeV  over  22.2  cm  to  give  an  accelerating 
gradient  of  39.2  MV/m.  With  this  result  and  the  observed  P6.3  MV/m  at  the 
gap,  we  see  that  for  this  case  a  transformer  ratio  R  =  9.85  has  been 
achieved . 

Several  interesting  aspects  of  this  simulation  should  be  noted. 

(1)  The  build  up  of  rf  in  the  cavity  is  of  a  transient  nature.  Were 
the  simulation  not  halted  at  t  =  30  ns,  the  field  amplitudes  would  increase 
beyond  the  observed  94.2  MV/m  until  limited  by  reflection  of  the  primary 
beam.  In  an  actual  device,  other  limitations  may  include  breakdown  in  the 
rf  structure,  losses  due  to  the  Q  of  the  cavity,  termination  of  the  primary 
beam,  or  acceleration  of  a  sufficiently  high  quantity  of  secondary  beam 
current . 
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(2)  The  conjectured  relationship  between  the  gap  field,  E 


and  the 


gap’ 

peak  axial  field,  E  .  ,  which  was  discussed  in  connection  with  Eq.  (2) 

v  axis  ‘  '  ' 

above,  does  not  hold.  Here.  v»  have  E  .  /'F.  =  1.67  and  J„(0)/J,,(kr.  )  = 

axis  gap  0  0  b 

2.51.  While  the  conjectured  relation  does  not  hold  in  a  precise  way,  it 

may  still  be  useful  as  a  qualitative  guide.  We  sill  expect  that  an 

increase  in  J()(0)/J^(kr^) ,  obtainable  by  increasing  r^,  will  result  in  an 

increased  E  .  /F.  .  This  will  be  investigated  below, 

axis  gap  & 

(3)  The  electric  field  of  56.3  MV/m  that  is  observed  across  the 
1.57  cm  gap  indicates  that  the  primary  beam  loses  0.883  MV  as  it  traverses 
the  gap.  This  energy  loss  is  verified  in  the  phase- space  plots  of  Fig.  8, 
where  the  primary  beam  particles,  which  have  0  <  z  <  10  cm,  are  deflected 
in  momentum  space  by  the  gap  voltage.  This  indicates  a  power  drain  of 
7.07  GW  at  1.27  GHz  and  is  sufficient  power  to  accelerate  secondary  beam 
current  in  the  500  A  range  over  this  short  distance.  With  a  longer 
accelerating  structure,  lower  currents  may  be  accelerated  to  higher- 
energies  . 

To  test  our  conjecture  that  higher  current  may  be  accelerated  to 

obtain  high  power  in  the  secondary  beam,  we  repeated  the  simulation  of 

Figs.  4-8  with  the  secondary  beam  current  increased  to  200  A.  We  found 

F,  -  51.6  MV/m  and  F.  .  =  91.3  MV/m.  Secondary  beam  particles,  injected 

gap  axis  r  J 

at  0.1  MeV,  were  accelerated  to  8.02  MeV  to  give  <E  .  >  =  35.7  MV/m  so 

axis 

that  R  -  9.78.  A  comparison  of  Jie.se  lesults  with  those  of  Figs.  4-8 
indicates  that  the  200  A  secondary  beam  does  not  significantly  load  the 
cavity.  We  also  see  that  1.58  GW  of  if  power  lias  been  transferred  from  the 
primary  to  the  secondary  beam. 

While  the  supposed  relationship  between  F.  and  F,  .  discussed  in 

gap  axis 

connection  with  Eq .  (2)  has  alieady  been  pi  oven  imprecise,  the  possibility 
of  obtaining  very  high  transformer  ratios  as  the  primary  beam  radius 
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approaches  r^  =  ^/k,  where  Jq(Jq  =  0,  remains  intriguing.  We 

investigated  this  by  repeating  the  simulation  of  Figs.  4-8  with  r^  =  8.0 

cm.  This  necessitated  an  increase  in  the  drift  tube  radius  to  r  =  8.4  cm, 

a  change  in  geometry  which  shifted  the  resonance  slightly  to  1.34  GHz.  At 

this  frequency,  ^/k  =  8.57  cm.  The  results  of  the  simulation  are  shown 

in  Figs.  9-12.  We  found  field  gradients  of  E  =  13.5  MV/m  and  E 

gap  axis 

34.4  MV/m.  Particle  plots  (not  shown)  indicated  that  the  secondary  beam 

particles,  injected  at  0.1  MeV,  were  accelerated  to  3.27  MeV  to  give 

<E  .  >  =  14.3  MV/m  so  that  R  has  been  increased  to  15.0.  Figures  9-12 
ax  1  s 

contain  the  following  results: 

(1)  With  r,  =  8.0,  we  have  E  .  /E  =  2.55,  an  increase  from  the 

'  b  axis  gap 

value  of  1.67  that  was  obtained  at  =  6.4  cm,  but  not  nearly  as  large  as 
J^COl/J^Ckr^)  =  11.9.  Note  that  the  transformer  ratio  was  similarly 
increased,  from  9.85  to  15.0.  As  stated  above,  we  have  only  a  qualitative 
ability  to  predict  results  as  r ^  is  changed. 

(2)  Figures  9  and  10  show  that  the  build-up  of  rf  fields  in  the  cavity 

is  of  a  transient  nature,  as  before,  but  much  lower  amplitudes  are  reached 

at  t  =  30  ns  than  in  the  r,  =  6.4  cm  case.  This  indicates  that  as  the 

b 

E  .  /E  ratio  is  increased,  the  interaction  between  the  primary  beam  and 
axis  gap 

the  rf  structure  is  weakened.  This  occurs  because,  at  a  higher  transformer 
ratio,  the  same  accelerating  field  in  the  rf  structure  gives  a  lower 
decelerating  field  at  the  gap  and  less  energy  is  drained  from  the  primary 
beam  per  cycle.  The  low  frequency  excitation  of  the  cavity,  apparent  in 
Figs-  4  and  9,  is  unchanged,  making  it  more  prominent  in  the  latter  case 
where  the  rf  fields  are  weaker. 

(3)  The  peak  electric  field  on-axis,  plotted  in  Fig.  10,  appears  to  be 
saturating  as  the  simulation  is  terminated.  It  is  not  clear  whether  this 
is  a  result  of  the  low  Q  of  the  numerical  structure  or  if  we  are  driving 
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the  cavity  slightly  off  resonance.  We  can  also  see,  from  Figs.  11  and  12, 
that  the  mode  structure  is  unchanged  from  the  previous  cases. 

Finally,  we  must  note  that  at  ry  =  8.4  cm,  f  =  1.34  GHz  is  very  close 
to  the  cutoff  frequency,  f£  =  ^c/2nrv  =  1.37  GHz.  In  a  practical 
device,  it  may  not  be  possible  to  increase  r^  and  ry  to  such  large  values 
at  this  frequency. 

IV.  Numerical  Effects 

To  understand  the  applicability  of  the  simulation  results  to  an  actual 

device,  it  is  necessary  to  examine  the  differences  between  such  a  device 

and  the  numerical  model.  Many  of  these  have  already  been  addressed.  One 

which  was  not  is  the  addition  of  a  center  conductor  to  the  drift  tube 

region  of  the  simulation  geometry,  which  provides  a  path  for  the  dc 

component  of  the  primary  beam  current.  The  significance  of  this  addition 

may  be  examined  by  considering  equivalent  circuit  elements  for  the  rf 

structure,  a  capacitive  load,  and  the  center  conductor,  an  inductive  load. 

These  elements  are  connected  in  parallel  and  are  driven  by  an  oscillatory 

current  source.  The  inductance  of  a  coaxial  line  varies  as  L  <=  log(rw/rc), 

where  r£  is  the  radius  of  the  center  conductor.  The  equivalent  circuit 

model  suggests  that  an  increase  in  rc  will  lower  the  inductive  load 

relative  to  the  capacitive  load,  lowering  the  voltage  across  the 

capacitance.  This  was  verified  by  increasing  the  radius  of  the  center 

conductor  to  r  =  5.0  cm  in  the  r,  =  6.4  cm  case.  This  had  the  effect  of 
c  b 

lowering  the  field  amplitudes  in  the  gap  and  on-axis  by  a  factor  of  1.7, 
but  left  the  transformer  ratio  unchanged.  Conversely,  the  circuit  model 
suggests  that  for  sufficiently  small  values  of  rc,  the  inductance  will  be 
so  high  that  it  will  behave  as  an  open  circuit.  In  this  ideal  case,  the 
entire  load  lies  across  the  capacitance. 
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To  discover  whether  or  not  the  radius  of  the  center  conductor  is 


sufficiently  small,  we  repeated  the  r^  =  6.4  cm  simulation  with  the  dc 

component  of  the  primary  beam  current  removed,  so  that  I^CkA)  = 

g( t)8sin(2nf t) ,  where  g(t)  is  an  envelope  function  as  before.  This  was 

accomplished  by  superimposing  an  appropriately  modulated  electron  beam  with 

a  dc  positron  beam,  and  allowed  us  to  compare  results  with  and  without  the 

presence  of  the  center  conductor.  With  the  center  conductor,  we  found 

E  =56.1  MV/m,  E  .  =92.6  MV/m  and  <E  .  >  =  38.8  MV/m,  comparable  to 
gap  axis  axis  Y 

the  results  of  Figs.  4-8.  Without  the  center  conductor,  we  found  E 

gap 

60.0  MV/m,  E  .  =  103.6  MV/m  and  <E  .  >  =  41.5  MV/m.  This  indicates  that 

axis  axis 

the  presence  of  the  center  conductor  reduces  the  fields  by  5-10  X. 


V.  Conclusions 

We  have  demonstrated  that  high  fields  and  transformer  ratios  can  be 
supported  by  a  MIREB-driven  accelerator,  with  several  interesting 
properties.  The  most  crucial  of  these  is  that  the  MTREB  is  so  strongly 
coupled  to  the  disk-loaded  rf  structure  that  power  in  excess  of  1  GW  may  be 
transferred  from  the  primary  to  the  secondary  beam,  despite  the  low  Q  of 
the  numerical  structure. 

We  have  found  that  the  build-up  of  the  rf  fields  in  the  structure  is 
transient  by  nature  and,  in  the  simulations,  peak,  accelerating  gradients 
were  limited  only  by  the  brevity  of  the  simulations.  In  an  actual  device, 
these  fields  will  continue  to  increase  in  amplitude  until  limited  by 
breakdown  in  the  rf  structure,  reflection  of  the  primary  beam  at  the  gap  or 
by  termination  of  the  primary  beam  pulse. 

We  have  also  considered  variations  of  the  geometry  to  successfully 
obtain  an  increased  transformer  ratio,  but  at  the  cost  of  weakening  the 
coupling  between  the  primary  beam  and  the  rf  structure.  We  have  also  found 
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that  the  conjectured  relationship  between  the  decelerating  field 
experienced  by  the  primary  beam  at  the  gap  and  the  peak  accelerating 
gradient  on-axis,  which  is  discussed  in  connection  with  Eq.  (2),  provides 
only  a  qualitative  guide  to  these  geometric  variations.  As  the  original 
conjecture,  contained  in  Ref.  2,  is  heuristic  in  nature  and  pertains  to  an 
idealized  physical  model,  this  is  not  a  surprising  result. 

Finally,  the  differences  between  these  simulations  and  a  practical 
experimental  configuration  have  been  discussed  in  some  detail,  suggesting 
that  similar  power  levels,  fields  and  transformer  ratios  may  be  obtainable 
experimentally. 
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Fig.  1.  MIREB-driven  accelerator  schematic. 
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Fig.  8  —  Particle  positions  in  phase-space,  y0,c  versus  z,  at  intervals 
left,  0<z  <  10  cm.  The  peak  energy  of  the  accelerating  secondary  beat 
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Fig.  8  (Continued)  —  Particle  positions  in  phase-space.  y0:c  versus  z,  at  intervals  of  0.2  ns.  The  primary 
beam  is  on  the  left,  0<<:  <  10  cm.  The  peak  energy  of  the  accelerating  secondary  beam  particles  is  noted  on 


Fig.  8  (Continued)  —  Particle  positions  in  phase-space,  yff.c  versus  z,  at  intervals  of  0.2  ns.  The  primary 
beam  is  on  the  left,  0<z  <  10  cm.  The  peak  energy  of  the  accelerating  secondary  beam  particles  is  noted  on 
each  plot. 
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Fig.  8  (Continued)  —  Particle  positions  in  phase-space,  y(3zc  versus  z,  at  intervals  of  0.2  ns.  The  primary 
beam  is  on  the  left,  0<z  <  10  cm.  The  peak  energy  of  the  accelerating  secondary  beam  particles  is  noted  on 
each  plot. 
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Fig.  11.  E  versus  z  on-axis  at  t  =  30.0  ns. 
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Fig.  12.  E  versus  r  near  the  point  of  peak  axial  field,  z  =  20.6  cm,  at 
z 


t  =  30  ns. 
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